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We investigate the orientational behaviour of a homeotropically ahgned nematic 
liquid crystal subjected to an oscillatory plane Poiseuille flow produced by an alternat- 
ing pressure gradient. For small pressure amplitudes the director oscillates within the 
flow plane around the initial homeotropic position, whereas for higher amplitudes a 
spatially homogeneous transition to out-of-plane director motion was observed for the 
first time. The orientational transition was found to be supercritical and the measured 
frequency dependence of the critical pressure amplitude in the range between 2 and 20 
Hz was in quantitative agreement with a recent theory. 

1 Introduction 

In nematic liquid crystals (NLCs) complex behaviour arises from the strong coupling between 
velocity v and director fi. For a steady flow along the x axis with a velocity field Vx{z), 
Vy = Vz = (rectilinear plane flow, typically of the Couette or Poiseuille type) the director 
will, in the absence of other torques, align in the flow plane {x — z plane) at the angle Oji = 



± arctg(y a3/a2) with the x axis if as/a2 > (the ± sign corresponds to positive/negative 
shear rate dv^/dz) Q]. In common nematics 0:3/02 is small but positive (~ 0.01), however, 
in some materials (in particular near a nematic - smectic transition) one has 03/02 < and 
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instead of flow alignment there is more complicated tumbling motion p, ^, 0. Note 
that a2 is always negative, so the sign of a-i/a2 is determined by that of as. 

When the velocity field oscillates periodically and symmetrically around zero [time av- 
erage {vx{z,t)) = 0] the situation becomes more complicated. Here we are concerned with 
oscillatory Poiseuille flow where the velocity field Vx{z^t) is induced by applying an al- 
ternating pressure difference AP cos(cjt) between the ends of a plane cell filled by a NLC 
(bounding plates parallel to the x — y plane). The planarly aligned case with the director 
oriented perpendicular to the flow plane, i.e., in the y direction, has been studied intensively 
in experiments by Pieranski and Guyon ^ and theoretically by Dubois- Violette and Man- 
neville [p!0| , [ll| . For very low frequencies a spatially homogeneous instability at some critical 
pressure amplitude has been found whereas for higher frequencies a transition to periodic 
rolls oriented along the flow direction was observed ||^, |^. 

When the director is prealigned within the flow plane the behaviour under oscillatory 
Poiseuille flow can be quite complex and has not yet been clarified completely. In the high- 
Ericksen number limit where elasticities (and therefore also boundaries) are neglected the 
time-averaged torque is zero and one has reversible behaviour (no attractors). Thus one 
expects the existence of a two-parameter family of solutions for the director oscillations de- 
pending on the initial orientation. When the effect of elasticities, which represent singular 
perturbations, is included and the shear rate is spatially nonhomogeneous (as is the case 
for Poiseuille, but not for Couette flow), there arizes an averaged torque, which produces 
attractors |]T2[. It is then found that, without surface anchoring, the flow- alignment solution 
is always unstable and one is left with one major attractor, which for 0:3/0:2 < is the orien- 
tation perpendicular to the flow plane, or, for 03/02 > 0, oscillations around an orientation 
close to it (in the y — z plane). In this case also the perpendicular orientation becomes 
unstable which is consistent with the results of [H, ^, |10|, |1T|. In both cases, the oscillations 
in the flow plane around the homeotropic direction are unstable. Including boundary effects 



one in fact expects a threshold for an out-of-pane transition ||T2 



In this letter we give the first experimental confirmation of this effect. We also present 
new, rigorous numerical calculations of the threshold as a function of frequency and compare 
with the experimental results. 



2 Experimental setup 

Experiments were carried out with the nematic MBBA [N-(4-methoxybenzylidene)-4-butyl- 
aniline], a flow aligning material (03, 02 < 0). The nematic was confined between two 3 
mm thick glass plates of lateral dimensions = 15 mm x Ly = 20 mm. The thickness 
d of the nematic layer was controlled by two polyester-foil spacers placed along the shorter 
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sides at y = ±Ly/2 (nominal thickness 23 fim). They also served to close the flow cell 
along these sides. The other sides were open. The nematic was held inside the cell by 
capillary forces. To induce homeotropic alignment of the director at the glass plates the inner 
surfaces were treated with the surfactant DMOAP [dimethyloctadecyl-(3-trimethoxysilyl)- 
propylammonium chloride] . 

The setup surrounding the cell was designed to apply an oscillatory air pressure difference 
Pi^2 = APi^2 cos(ci;t + (pi^2) to the two open sides of the cell (fig. |I]). The cell was mounted 
inside a brass block of outer dimensions 4 x 5 x 10 cm. The open sides of the cell were 
in contact with channels in the brass block which widened towards the cell in order to 
gain the width Ly. As a result the pressure was distributed rather evenly over the entire 
width of the cell. The channels connected, via cavities of adjustable volume, to two vertical 
cylinders accommodating pistons driven sinusoidally in opposing sense by an electromotor. 
The displacement was f ^ = 5 cm'^ and the frequency / = uj/{27c) ranged from 1 to 20 Hz. 
The volume of the cavities could be adjusted by screws in order to regulate the "passive" 
volume Vi of air between the lower position of the pistons and the cell. This allowed to vary 
the amplitude AP of the applied pressure in a range from 1 to 20 kPa. Pressure sensors were 
connected to the open sides of the cell in order to measure the pressure head with respect 
to the surrounding air. 

The response of the nematic was detected optically by exploiting its birefringence. We 
used a light source with maximum intensity at wavelength 610 nm (cold source with filter) 
and crossed polars (polariser || e^, analyser || iy). The transmitted light intensity was 
measured by a semiconductor detector, either integrating over an area of 28.3 mm^ or with 
a spatial resolution of 2 fim (long-range microscope). The photo detector in the microscope 
could be replaced by a CCD camera. Interfaces for the signal processing were either an A/D 
board or a frame grabber with a resolution of 512 x 512 pixels and 256 grey scales. The 
temperature of the brass block was kept at 25.0 ± 0.1°C by a water thermostat. 

Due to the large aspect ratio of the fiow cell {Ly/d ^ 10'^) one may expect rectilinear 
Poiseuille fiow {vy = = 0) except in thin boundary layers near the edges. In order to test 
the fiow field and to calibrate the pressure sensors the motion of AI2O3 tracers (diameter 
5 fim) was recorded with the CCD camera and analysed. No y component of the fiow was 
detected. 

An alternating pressure applied to a fiuid layer decays by compressibility effects over a 
length 6 = d\j2pc^ / {n'^ur]) |l^, |T^. Here p is the mass density and c is the speed of sound. 
The viscosity for the homeotropic orientation is given hj rj = (— ^2 + «4 + «5)/2. Using 
MBBA material parameters (p = 10^ kg/m^, c = 1.5-10^ m/s [0, rj = 135.4-10"^ N-s/m^), 
cell thickness d = 23 pm, and the frequency / = 20 Hz one obtains 5 ~ 12 cm which 
is considerably larger than the cell length = 15 mm. Moreover, since in our setup the 
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pressure difference is applied symmetrically at both sides of the capillary the flow is expected 
to be practically homogeneous along the x direction. This was also verified by tracing the 
AI2O3 particles. 



3 Experimental results and comparison with theory 

Let us define the director components as n = (sin ^ cos 0, sin 6' sin 0, cos6') with the polar 
angle 6 measured with respect to the z axis and azimuthal angle measured from the x 
axis. For amplitudes AP below any instability the director is expected to oscillate within 
the flow plane, i.e. = 0, and in our geometry of crossed polars the intensity should be 
minimal (and constant in time). As a result of an out-of-plane transition the intensity should 
increase except when = 7r/2 and = vr. In fig. |^ we show examples of the signals of the 
two pressure sensors together with the signals of the transmitted light intensity measured 
by the photo diode as functions of time for an amplitude below and above threshold. We 
have 1^2 = 1^1 + IT and APi = AP2 = AP/2. 

The signal of the light intensity in fig. ^ exhibits two rather sharp double peaks over 
the oscillatory flow period. We interpret the minimum as the (rapid) passage of the director 
through = 71 /2. The long low-intensity interval between the peaks then results from 
approaching and vr while the director reverses its motion (trajectories of the director motion 



are shown in ref. |jT2[). Slightly above the out-of-plane instability threshold one obtains for 



the azimuthal angle the following expression |T2[ 

tan0 = x/[8a(z/(i) sintut], (1) 

where x ("^ 1) is the out-of-plane angle between the director and flow plane {x—z plane) when 
the director passes through = 7r/2 and a = A/d is the dimensionless flow displacement 
amplitude at midplane. Clearly |0| is maximal at the cell boundaries z = ±d/2 and the 
leading contribution to the transmitted light intensity comes from there (for x <^ 1). This 
gives an approximate time dependence of the transmitted light intensity 



iTr 



sin2[20(;z = d/2)] 



X4:a sin oot 



(2) 



+ {Aa sin Luty 

Equation @ describes the observed signal of the light intensity (fig. ^) and agrees with the 
results of full simulations of the optical response based on the numerical solution of Maxwell's 
equations with the director distribution calculated from the nematodynamic equations (see 
below). In principle such type of signal could also be caused by an elliptical motion of the 
director without breaking the left-right symmetry. We have excluded this possibility by 
measurements with the analyser rotating with the driving frequency. 
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Since the out-of-plane transition involves the breaking of a two-fold symmetry one expects 
the appearance of patches of the two states. Indeed, under the microscope we always observed 
domain boundaries indicating patches of ~ 0.1 — 1 mm^ area (fig. After some initial 
coarsening they did not change much over times of several minutes, except for oscillatory 
motion of the boundaries with the flow field (see the bottom part of fig. 

In order to obtain quantitative results for the threshold we determined the average in- 
tensity of the transmitted light {Irrit)). In fig. ||a {Ixrit)) is shown versus the pressure 
amplitude for a frequency / = 10 Hz. By integrating eq.(^ one finds {Ixrit)) ~ x- Since 
one expects x ~ V ~ "we have fitted the measured values of {Itt) accordingly. The 



best fit curve {Itt) = 16.4a/AP — 7.3 is shown in fig. ^a (solid). Clearly the out-of-plane 
motion sets in at a well-defined threshold APc via a continuous (forward) bifurcation (no 
hysteresis), as predicted by theory 



Direct numerical simulations of the standard set of nematodynamic equations [T^ where 
the director n and velocity v are functions only of the distance z from the boundaries and 
time t (see ||13[ for details) were performed using finite differences for the spatial derivatives 
and the predictor-corrector scheme for the time discretisation. The critical amplitude of the 
pressure gradient APc(/)/-^x was calculated for the frequency range 0.5 Hz < / < 20 Hz, 
layer thickness c? = 23 /xm and MBBA material parameters at 25°C |1T8| , [T9[| (elastic constants 
in units of 10~^^ N: Ku = 6.66, K22 = 4.2, and K^a = 8.61; viscosity coefficients in units 
of 10-3 N-s/m^: ai = -18.1, as = -110.4, 03 = -1.1, = 82.6, 05 = 77.9, ag = -33.6; 
mass density p = 10^ kg/m^). The results of the numerical computations together with the 
measured values of APc(/) are shown in fig. The agreement of the measured values with 
the calculated ones is quite good (no adjustable parameters). Discrepances near the ends of 
the frequency interval can be explained by deviations from the sinusoidal driving in these 
regions. One notices that for not too low frequencies {oot^ ^ 1, where is the director 
relaxation time) APc rises linearly with u. This corresponds to the critical Poiseuille flow 
amplitude Oc being independent of u, which can be understood from simple scaling arguments 
fl^ . Quantitative discrepancies with our earlier theoretical work results from the fact that 
there the effect of the director distribution on the flow field, which leads to deviations from 
the parabolic profile, was not included. 



4 Conclusion 

For the first time the out-of-plane orientational transition in homeotropically oriented ne- 
matic liquid crystal under oscillatory Poiseuille flow was observed experimentally. The phe- 
nomenon is interesting because the underlying torque on the director resulting from spatially 
inhomogeneous driving in the presence of diffusive coupling, arises quite unexpectedly. 
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The analysis of the transmitted hght intensity shows that the instabihty has a supercriti- 
cal character without hysteresis. The frequency dependence of the critical pressure amplitude 
is in a good quantitative agreement with theoretical results without any adjustable param- 
eters. Further efforts will be made to increase the range of frequency and pressure of the 
measurements. The theory predicts that at high driving amplitudes the director saturates 
at an orientation perpendicular to the flow plane (0:3 > 0) or near to it (03 > 0) unless 
a further transition to a spatially periodic state intervenes. 

Let us finally mention that a more complicated scenario is expected for boundary an- 
choring corresponding to the flow-alignment angle or, alternatively, simple planar anchoring 
of the director with an additional magnetic field in the flow plane at an angle of 45° with 
respect to the x axis. Then a transition to a slow time-periodic director precession has been 
predicted theoretically |]T^. This director motion emerges with increasing flow amplitude 
through a homoclinic bifurcation and disappears through a Hopf bifurcation. It would be 
interesting to perform experiments to test this prediction. 
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electro motor 




Figure 1: Cross section of the experimental setup. 



8 




Figure 2: Signal of the pressure sensors (dashed and dot-dashed hues) and transmitted hght 
intensity (sohd hne) for a subcritical pressure (a) and a supercritical pressure (b). 
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Figure 3: Snapshot of patches with left- and right-side out-of-plane director configurations 
(upper part) and x — t plot showing time oscillations of domain boundaries (lower part) far 
from threshold. The crossed polars were rotated slightly away from the symmetric orientation 
to obtain different intensities in the two configurations. The dark vertical stripes pertain to 
dust particles on the glass plates. 
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Figure 4: a) Transmitted light intensity averaged over the flow period for increasing (O) and 
decreasing (A) pressure amphtude (/ = 10 Hz), b) Frequency dependence of the critical 
pressure amplitude APg for out-of-plane transition: experimental (squares) and calculated 
(solid line). 
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